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Far too often are multiparticle final states studied and models tested on merely 
single-particle spectra and their integrals, the average multiplicities: A multipar- 
ticle final state is a non-linear, complex system and the essential information on 
the dynamics is contained in the particle correlations and their integrals measuring 
the fluctuations! 



1 The formalism 

We start by defining symmetrized inclusive g-particle distributions 

q 

p q (p 1 ,...,p q )=a^ lda q (pi,...,p q )/Y[dp q , (1) 

l 

where cr q (j>i, ■ ■ ■ ,p q ) is the inclusive cross section for q particles to be at 
Pi, . . . ,p q , irrespective of the presence and location of any further particles, pi 
is the (four-) momentum of particle i and Otot is the total hadronic cross sec- 
tion of the collision under study. For the case of identical particles, integration 
over an interval Q in p-space yields 



dpi . . . / dpgPgipx, ...,p q ) = (n(n - 1) . . . (n - q + 1)) , (2) 
in isi 

where n is the multiplicity of identical particles within f2 in a given event and 
the angular brackets imply the average over the event ensemble. 

Besides the interparticle correlations, the p q in general contain "trivial" 
contributions from lower-order densities. We, therefore, consider a sequence 
of cumulant coorelation functions C q (pi, . . . ,p q ) which vanish whenever one 
of their arguments becomes statistically independent of the othersu 

C 2 (Pl,P2) = P2(Pl,P2) ~ Pl{Pl)Pl{P2) , 

C 3 (pi,P2,P3) = P3(P1,P2,P3) -^2pi(l)p 2 (p2,P3) + 2p 1 (pi)pi(2p ) p 1 (p 3 ) , 

(3) 

etc. Deviations of these functions from zero shall be addressed as genuine 
correlations. For recent discussions, it is important to note that not the p q 
but the C q are additive in case of N independent overlapping samples. 
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It is often convenient to normalize the functions p q and C q by the uncor- 
rected densities, 

Rqipi, . . ■ ,Pq) = Pq(p q , . . . ,Pq)/pi(pi) ■ ■ -Pl(Pq), (3) 
K q (pi, . . . ,p q ) = C g (pu . . . ,Pq)/pi(pi) ■ ■ -Pl{Pq)- (4) 

Important relations between the two sets of functions are: 

R 2 = 1 + K 2 ; R 3 = 1 + 3K 2 + K 3 (5) 
Because of the additivity of C q , the K q will show a dilution according to 
Kq N ^ = N( 1 ~ q ) Kq 1 ^ for N independent overlapping sources. 

Correlations have been studied extensively for q — 2 and a strong positive 
correlation is found for p\ « p 2 n To be able to do that for q > 3, one must, 
however, work via the integrals Eq.(2), in limited phase-space cells. B 

In practical work, with limited statistics, it is almost always necessary 
to perform averages over more than a single phase-space cell. Let f2 m be 
such a cell (e.g. a single rapidity interval of size 5y) and divide the phase- 
space volume fi into M non-overlapping equal cells fl m of size SQ = Q/M. 
Normalized cell-averaged factorial moments 13 are then defined as 



M J&a, Pq(P\i ■ ■ ■ iPq) IILl d P' _ 1 V~> ("mKn-1) ■ ■ ■ (n m -q+l)) 



A/ i 



(Ssy P^ d P 



1 l.s..PU>)up) ^ 

r (6) 

Likewise, cell-averaged normalized factorial cumulant moments K q (Sy) 
may be defined □ by replacing the p q by C g in (g). In analogy to Eq. (||), they 
are related to the factorial moments by 

^2 = 1 + ^2 ; F 3 = 1 + 3K 2 + K 3 . (7) 

Essential porperties are: 

1. The Poisson-noise suppression of the F q , contrary to e.g. (n q ) /(n) q . 

2. Independent emission leads to F q = 1 and K q = for all q. 

3. The factorial moments resolve the high-n m tail of the multiplicity dis- 
tribution and are, therefore, sensitive to high-density fluctuations. 

A fruitful further development □ is the use of integrals of p q or C q over 
a strip domain, rather than a sum of the box domains. One not only avoids 
unwanted side-effects such as splitting of density spikes, but also drastically 
increases the statistical significance at a given resolution. It, furthermore, 
allows to work in terms of inter-particle distance measures, as the (invariant) 
four-momentum difference Q = y— (Pi — Pi) 2 ■ 

As shown in 13, a "smooth" (rapidity) distribution, which does not show 
any fluctuations except for the statistical ones, has the property of F q (Sy) 
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Figure 1. Comparison of density integrals for qn= 2 as a function of 2 log(l/Q 2 ) for e+e _ 
(DELPHI) and hadron-hadron collisions (UAl).Q 

being independent of the resolution Sy in the limit Sy — > 0. On the other 
hand, if self-similar dynamical fluctuations exist, the F q obey the power law 

F q (5y) oc , (<5y->0). (8) 

The powers g (slopes in a double-log plot) are related i to the anomalous 
dimensions d q — 4> q /(q — 1), a measure for the deviation from an integer di- 
mension. Equation (0) is a scaling law since the ratio of the factorial moments 
at resolutions L and £, R — F q (£) / F q (L) — (L/£)^ q , only depends on the ratio 
L/£, but not on L and £, themselves. 

The scaling behavior is related to the physics of fractals. Fractal behavior 
is indeed expected for multiparticle production both from QCD branching and 
from a QCD phase transition. 

2 The State of the Art on Power-Law Scaling (Intermittency) 

2.1 Existence of scaling: The suggestion that the F q or K q might show power- 
law behavior has spurred a vigorous experimental search and (more or less! 
linear dependence of In F q on — In SQ has been found in all types of collisions B 
Recent e + e~ results will be presented here by De Wolf. 

Of particular interest is a comparison of hadron-hadron to e + e _ results in 
terms of same and opposite charges of the paxticles involved. This comparison 
has been done for UA1 and DELPHI data in 13 and is shown in Fig. 1 for q = 2. 
Important differences between UA1 and DELPHI can be observed: 

For relatively large Q 2 (> 0.03 GeV 2 ), where Bose-Einstcin effects do 
not play a major role, the e + e _ data increase much faster with increasing 
2log(l/Q 2 ) than the hadron-hadron data. The increase in this Q 2 region is 
very similar for same and for opposite charges. At small Q 2 , however, the 



SAME CHARGES x,l' 



} DELPHI ■ UA! 
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e + e~ results approach the hadron-hadron results and the increase is faster for 
same than for opposite charges. So, at least two processes are responsible for 
the power-law behavior: Bose-Einstein correlation at small Q 2 following the 
evolution of jets at larger Q 2 . 

As such, Bose-Einstein correlations from a static source are not power- 
behaved. A power law is obtained a if i) the size of the interaction region is 
allowed to fluctuate, and/or ii) the interaction region itself is a self-similar 
object. 

2.2 Multifractal or Monofractal: Fractals with all d q equal are called monofrac- 
tals, those with d q depending on q are called multifractals. A second-order 
phase transition (e.g. from a QGP) would give rize to a monofratal behavior 
of particle production, a cascade process (as e.g. QCD branching) leads to a 
multifractal behavior. Even though the d q and their q dependence are found 
to be largest in e + e~ and smallest in heavy-ion collisions, so far all types of 
collisions (including AA collisions) exhibit multifractal behavior. We shall 
hear more about monofractal behavior from N. Antoniou. 

2.3 Self-similar or Self-affine Fractal: If the power law (|8|) holds when p- 
space is partitioned by the same factor in different directions, the fractal is 
called self-similar. If, on the other hand, (^) holds and only holds when space 
is partitioned by diffcrcnt_Eactors in different directions, the corresponding 



Self- affinity can be characterized by the Hurst exponents Jid Hij = 
laMi/lnMj (0 < < 1) with M % (i = 1,2,3; M x < M 2 < M 3 ) be- 
ing the partition numbers in the phase-space variables pi. The Hij can be 
obtained from the saturation of the one-dimensional hiF2(Spi) distributions. 

For hh collisions, the Hij for kmgitudinal-tpjisverse combinations is in- 
deed smaller than unity for NA22 E3 and NA27 £3, while it is consistent with 
unity within the transverse plane (see Wang Shaoshun for hh and Chen Chang 
for first results on e + e~ collisions) 

2-4 ■ QCD + LPHD have been quite successful in the description of single- 
particle £ distributions and the energy evolution of its maximum position £* 
and integral n. According to Eqs. ([!]) and (||), these correspond to q = 1 only 
and, therefore, neglect any vital information on particle correlations, 
a. High-Order Correlations: The ratio H q = K q /F q reflects the genuine 
g-particle correlation relative to the g-particle xlfinsity. They have been cal- 
culated for gluons at different orders of QCD.E3 For NNLLA, H q decreases 
to a negative first minimum for q = 5 and shows quasi-oscillation around 
for q > 5. Assuming the validity of LPHD, such a behavior is also expected 
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for the charged particles. For e + c~ — > hadrons at y/s — Mz, SLD LLj indeed 
observed such a behavior. Howevet, the same behaviour is observed iip-hh 
collisions between 20 and 900 GeV O and even in hA and A A collisions S3 

An L3 analysis of sub-jet multiplicities tB reveals that this behavior ap- 
pears only for energy scales < 100 MeV, far away from the perturbative region. 
Furthermore, similar oscillatory behavior is observed for a large variety of MC 
models, including independent fragmentation. One is tempted to conclude, 
that at present energies, the oscillations of H q are unrelated to the behavior 
predicted by the NNLLA calculations. They may be related to the energy 
conservation first approximated in NNLLA, but what do we have to think 
of a theory approximating energy conservation only in such a high order of 
perturbative expansion? 

The oscillations are reproduced by a superposition of 3 negative binomials 
corresponding to two-jet, mercedes-like three-jet and intermediate events, 
respectively.Ea So, the oscillations originate from an interplay between hard 
(jet topology) and soft (fragmentation) phenomena. 

b) Local fluctuations: In perturbative QCD, the fractal structure of jets follows 
from parton branching and the <j> q are directly related u to the anomalous 
dimension 70 = (Gcts/ir) 1 / 2 , at lowest order via d q = D — j (q + l)/q, where 
D is the topological dimension of the analysis and the second term corresponds 
to the Renyi dimensions D q of the fractal. 

Analytical QCD predictions have indeed been given for fluctuations in 
angular intervals. a Assuming LPHD onaiagaiiL^these predictions can be com- 
pared to experimental data. DELPHI B, L3 E3 and ZEUS B show that the 
predictions disagree with the data. Energy-conservation effects are indeed 
sizeable.E.3 Recent ZEUS results will be presented here by L. Zawiejski and 
the state of understanding will be summarized by W. Ochs. 

2.5 Event-to- Event Fluctuations, Erraticity, Voids and Entropy: The F q as 
defined in (^|) contain averages over all bins within one event and over the 
whole event sample. Event-to-event fluctuations lost in this procedure can 
be quantified c3 by defining the F q per event and determining moments of 
the distribution P{F q ) as C q , p — ((F q ^l$ / (F q ) p . However, in an analysis of 
the influence of statistical fluctuations JH3 it turned out that these completely 
dominate, even in heavy-ion collision. 

When the event multiplicity is low, the gaps between neighboring particles 
carry more information than multiplicity spikes. The—method has therefore 
been extended to measure the rapidity gaps (voids). E3 Experimentally the 
erraticity measures have been studied on NA27 B3 (see Wang Shaoshun). 

A method most sensitive to the maximum of the probability distribution 
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is a measure of entropy suggested in Ea. 

Furthermore, results based on event-to-event pt fluctuations will be pre- 
sented here by Bai Yuting and T. Sugitate. 



3 Bose-Einstein Correlations 

3.1 Alternative Views: As can be seen on Fig. 1, BEC play a role at small 
Q, both in hh and e + e~ [ _colhsions. Conventionally, they are studied within 
intensity interfcrometryH3 Basic ingredients are incoherent (chaotic) boson 
emission from a spherically symmetric static source and production ampli- 
tudes independent of momentum. Undere these conditions, 

R 2 = l + K 2 = l + \p(Q)\ 2 (9) 
R 3 = 1 + |pi 2 | 2 + |p 23 | 2 + |p 3 i| 2 + 2Re(pi 2 p 23 p3i) , (10) 

where p(Q) is the Fourier transform of the space-time density of boson emit- 
ters. L. Gurvits will show us, however, how far astrophysics has developed 
away from intensity to amplitude interferometry, in the meantime. 

A formalism particularly handy for the fully-dimensional_treatment of a 
dynamical emitter is the so-called Wigner-function formalismtJ This is based 
on the emission function S(x,p), a covariant Wigner-transform of the source 
density matrix. S(x,p) can be interpreted as a quantum-mechanical analogue 
of the classical probability that a boson is produced at a given space-time 
point x — (t, r) with a given momentum-energy p — (_E, p) r^ee T. Csorgo). 

Alternatively, BEC are introduced into string models S3 An ordering in 
space-time exists for the hadron momenta within a string. Bosons close in 
phase space are nearby in space-time and the length scale measured by BEC 
is not the length of the string, but the distance in boson-production points 
for which the momentum distributions still overlap (B. Andersson). 

3.2 Recent experimental results: Bose-Einstein correlations are by now a well 
established effect in the hadronic final states of Z and W decay (intra-W 
BEC).cJ The important question is that of BEC between pions each originat- 
ing from a different W in fully hadronic W + W~ final states (inter-W BEC). 
If existent, such a correlation could, on the one hand, cause a bias in the mass 
determination of the W. On the other, it could serve as a pion-interfcrometry 
laboratory for the measurement of the space-time development. The recent 
status of the search for inter-W BEC will be covered by J. van Dalen. 

Detailed knowledge of intra-W BEC would be necessary to understand 
inter-W BEC. Limited statistics prevents that, but a large amount of high 
statistics information has recently been obtained for the Z and other reactions. 
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For shortness, I am forced to refer to recent Zakopane lectures La on this topic 
and to limit myself to just listing the main results: 

Bose-Einstein correlations exist also in e + e~ collisions and DIS, with ra- 
dius parameters independent of y/s and Q 2 , respectively (E. De Wolf and L. 
Zawiejski). 

The region of homogeneity is elongated along the event axisEl This is not 
reproduced by JETSET + LUBOEI, but expected fromB. The short range 
of this region is evidence for a strong momentum-position correlation. 

The correlator Ki is steeper than Gaussian. EE This is contrary to most 
of the-assumptions and interpretations in the literature, but, again expected 
from H. 

All radius parameters show a l/y/m T dependence in all types 
reactions.EH This may be explainable already from the Heisenberg principle 
but is expected from an inside-outside cascade with transverse flow £3 More 
results from RHIC here from T. Sugitate. 

In hh and AA collisions, there is evidence for dilution of K 2 due to in- 
dependence of overlapping mechanisms .u£] Search for a similar dilution due to 
overlap of string pieces in gluon jets will be presented by N. van Remortel. 

The pion ^mission function in space-time has been extracted for hh and 
AA collisionsO In the transverse plane, it is Gaussian for AA collisions but 
ring-shaped (!) for hh collisions. Study of the Z is underway (T. Csorgo). 

For 7r°7r°, according to the string model, the radius should be smaller 
than for charged pions. Wes Metzger will present the experimental evidence. 

Genuine three-particle correlations contain additional information on a 
phase or, alternatively, on dilution of K q proportional to A^ 1 " 9 ) for N inde- 
pendent sources. W. Metzger will give recent results on the Z. 

Modelling these observations in a Monte Carlo code is a major challenge 
(see G. Wilk, K. Fialkowski). 

4 Conclusion 

We need a model containing parton branching and energy conservation plus 
soft fragmentation to reproduce the fractal behavior of particle production. 
JETSET and ARIADNE are doing surprisingly well there, HERWIG a bit less. 
We further need the model to do away with old prejudice from conventional 
pion interferometry. A large amount of new information is (becoming) avail- 
able on BEC at the Z. Models neglecting this information I would consider a 
loss of time. 
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